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Abstract 

We address the question: Given an arbitrary initial state and a general physical interaction what 
is the minimum time for reaching a target entangled state? We show that the minimum time is 
inversely proportional to the quantum mechanical uncertainty in the non-local Hamiltonian. We 
find that the presence of initial entanglement helps to minimize the waiting time. Furthermore, we 
find that in a bi-local rotating frame the entangling capability is actually a geometric quantity. We 
give an universal bound for the time average of entanglement rate for general quantum systems. 
The time average of entanglement rate does not depend on the particular Hamiltonian, rather on 
the fluctuation in the Hamiltonian. There can be infinite number of nonlocal Hamiltonians which 
may give same average entanglement rate. We also prove a composition law for minimum time 
when the system evolves under a composite Hamiltonian. 
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I. INTRODUCTION 



Quantum entanglement plays a pivotal role in the emerging field of quantum information 
theory [l[ . Creation, storage and processing of entangled states are challenging experimental 
tasks in any quantum information processing devices 0, 0]- In recent years there have 
been considerable progress in understanding various types of entangled states. However, in 
multiparticle situations, many questions are still unexplored. We know that if we start from 
an unentangled state of two (or more) particles and allow them to interact, then the state 
may evolve to an entangled state depending on the initial state and the type of interaction 
Hamiltonian. The Hamiltonian that is capable of creating entanglement is not sum of 
local Hamiltonians, rather it involves a nonlocal part. This is so, because sum of the local 
Hamiltonians will always give rise to tensor product of local unitary evolution operators and 
we know that the later cannot create any entanglement. In some situation, one would like 
to generate a particular entangled state for some specific quantum information processing 
task. For example, if one could create a maximally entangled state in a controlled way, then 
that can be used for quantum teleportation, dense coding, remote state preparation and so 
on. However, what type of initial state and what kind of non-local Hamiltonian can give 
the desired entangled state is far from clear. Also, it is not the case that immediately after 
switching on the interaction Hamiltonian we will have the target entangled state. 

There has been considerable interest in the study of dynamics of entanglement. In par- 
ticular, a pertinent question in this context is given a nonlocal Hamiltonian, how well one 
can create entangled states. This has been address by Diir et a/ (sf. They have introduced 
entangling capability of Hamiltonian for two qubits and also generalized for two qudits. It 
was found that it is better to start with initial entangled state and the best initial entangle- 
ment is independent of the physical process. Also, it was shown that one can improve the 
capability if we allow fast local operations and ancillas (in some cases) [sj. 

In this paper we address the following question: Given an arbitrary product state and 
a general physical interaction what is the minimum time for creating a target entangled 
state? This is precisely the entangled brachistochrone problem. This question is not only of 
fundamental interest but also of practical importance. The answer will be relevant in many 
experimental context where we want a fixed entangled state after having the interaction 
on. Because, that will allow us to generate the desired entangled state in a controlled way. 
This will also answer the question what is the minimum time an experimentalist needs to 
wait to create the target entangled state. In the course of our investigation, we show that 
the entanglement capability of non-local Hamiltonains is a geometric quantity. We prove 
a universal bound on the time average of the entanglement rate. Our result suggests that 
even though we cannot tell whether entangling rate and disentangling rates are same, but 
the maximum value of the time average of entanglement and disentanglement rates are the 
same. 

We prove that the minimum time for creating a fixed entangled state not only depends 
on the initial state and final state, but also depends on the speed at which system evolves. 
The speed of the quantum evolution is governed by the fluctuation in the Hamiltonian Aif. 
In particular, we show that the minimum time required to reach a target entangled state 
depends inversely on the quantum mechanical uncertainty in the non-local Hamiltonian. We 
illustrate this for two qubits. Further, we will show that the time average of entanglement 
rate is upper bounded by 21ogA^^^, where N is the number of Schmidt coefficients in the 
entangled state, AH is the fluctuation in the nonlocal Hamiltonian and 5*0 is the shortest 
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geodesic connecting the initial and the final states. As an application of this bound we 
estimate what could be the ultimate entanglement rate for non-local Hamiltonians. 

The present paper is organized as follows. In Sec. II, we briefiy discuss the geometric 
uncertainty relation and the bound for the minimum time for entangled states. In Sec. 
Ill, we illustrate the bound for two qubits. We find that initial entanglement can help to 
reduce the waiting time in reaching a maximal entangled state. In Sec. IV, we show that 
entanglement rate in a bilocal rotating frame is actually related to the quantum fiuctuation 
in the nonlocal Hamiltonian (the speed of the quantum evolution). Also, we provide a 
universal bound for the time average of the entanglement rate. In Sec. V, we discuss the 
question of minimum time when the system is driven by composite Hamiltonian. Finally, 
we conclude our paper in Sec. VI. 



II. GEOMETRIC UNCERTAINTY RELATION AND MINIMUM TIME 

Here, we consider bipartite systems but some of the results will be valid for more general 
situations. Whenever the result holds only for bipartite systems, we will say so. Consider an 
arbitrary initial product state |^I^(0)) = |V'(0)) ® |0(O)) of a bipartite composite system. The 
nonlocal Hamiltonian for the system is given by H = I + 1 ^ H2 + H^^t. Under the action 
of the Hamiltonian the state evolves unitarily as |^I^(0)) |^(^)) = exp(— zift/^)|\l/(0)). 
Now, depending on the type of initial state and Hi^t the resulting state will be entangled. 
Suppose, we evolve the system from |^(0)) |^(^)) and our target entangled state is 
I^^t) = |^(^))- Then, the question is what is the minimum time Tmin for which we should 
evolve the combined system to get the desired state? 

To answer this, we need some geometric ideas. Let be a set of vectors in Ti. = 7Yi®7-^2- 
If these vectors are not normalized we can consider a set of vectors |\I'| |} of norm one 
in C. The set of rays of H is called the projective Hilbert space V(Hi ®'H2)- If dim?-^i = 
and dim7^2 = M, then H ~ C^^^ The projective Hilbert space is P = (C^*^ - {0})/[/(l) 
which is a complex manifold of dimension {NM — 1). This can also be considered as a real 
manifold of dimension 2{NM — 1). Any quantum state at a given instant of time can be 
represented as a point in V via the projection map H : — >■ |^')(\&|. The evolution of the 
state vector can be represented by a curve C : t ^ ^ whose projection H(C) lies 

in V. Here, smooth mappings C : [0,t] — > £ of an interval into a differentiable manifold are 
called smooth curves in the given manifold f^. 

Now consider the unitary time evolution of a bipartite quantum system |^(0)) —>■ = 
exp{—iHt/h)\'${0)), where H is the nonlocal Hamiltonian. If the state at later time is 
entangled, then using the Schmidt decomposition theorem we can write the combined state 
as 

N 

\^{t)) = J2^/Kj^\anmbn{t)), (1) 

n=l 

where A„(t)'s are the Schmidt coefficients with ^„ A„(t) = 1 for all time; \an{t)) and \bn{t)) 
are the orthonormal Schmidt basis. As the state evolves, it traces a path in the state space 
V. There is a natural metric on V which is the Fubini-Study metric, defined as [g], 0, [q] 

A H"^ 

dS' = A{l-\{^{t)\^{t + dt))\')=4^dt\ (2) 
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where AH^ = (\I/(t)|ff^|\E'(t)) — is the quantum mechanical uncertainty in 

the nonlocal Hamiltonian. For time-independent Hamiltonians AH is also time independent 
(because H and U{t) commutes). Using the notion of infinitesimal distance, one can define 
the rate or the speed V at which the quantum system evolves in the projective Hilbert space. 
This is given hj V — ^ — We can write down an explicit expression for the speed as 

(with the initial state as the product state) 

= ^[AH', + AHi + AHl, + C{H^, H,^,) + C{H,^t, H^) + C{H^, H,^,) + C{H,^,, H^MS) 

where AHf = {miHfim) " {m\Hi\m)', = imiHiim) - immmr, 

and C{A,B) = (^(0)|A5|^(0)) - (^(0)|A|^(0))(^(0)|S|*(0)) being the correlation be- 
tween A and B (here A and B can be Hi, H2 or i/jnt)- 

We can also give an expression for the speed of the entangled state evolution using the 
Schmidt basis and coefficients. Now, using (1) the Fubini-Study entangled metric can be 
expressed as 

dS^ = 5]A„(i)((a„(i)|d„(^)) + (6„(i)|6„(i)))di2 

n 

- [5^A„(t)(z(a„(t)|d„(t))+i(6,(t)|6„(t)))]' df 

n 

- 2^ VAn(t)A^(t)(a^(t)|d„(t))(6„(t)|6n(t)) + E (4) 

This metric is U{1) gauge invariant and is independent of the detailed dynamics of the 
systems. This is so, because, there can be many non-local Hamiltonians which may give the 
same path in V. If the state is not entangled (one of the Schmidt number is one and others 
are zero), then |^'(t)) may be written as |^(t)) = \a{t))\b{t)) and the metric is then given 
by dS"^ = dSf + (iS'|, where dSf, {i = 1,2) are the Fubini-Study metrics for the individual 
subsystems. In the above expression, the second and third terms represent the effect of 
entanglement on the metric structure. 

Let us apply the nonlocal Hamiltonian for a time period T and consider the quantum 
evolution |^(0)) |^(T)). Suppose that the state |^(T)) is our target state. The total 
distance travelled by the quantum system in going from the initial state to the target state 
(as measured by the Fubini-Study metric) is given by 5" = 2 J^, ^ dt — 2Vt, where V is 
given in (3). This is typically a longer path in the projective Hilbert space of the entangled 
system. Now, consider the shortest path 5*0 between the initial and the target state. It is 
given by 

\{^m^{T))\'^cos'^. (5) 

Since the actual distance S is greater than or equal to the shortest distance 5*0 (the Anandan- 
Aharonov version of the uncertainty relation) we have 

AH T>^. (6) 
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From the above geometric relation we have a bound on the time required to reach a target 
entangled state which is given by 



T>^cos-'\{^m^im. (7) 

Therefore, to minimize the waiting time, we have to evolve the system faster and minimise 
the shortest path. However, if the initial and target states are fixed, then the waiting time 
can be reduced by evolving the system faster. It may be noted that the above equation is 
valid both for qubits and qudits as well as for multiparticle quantum systems. 



III. OPTIMAL TIME FOR TWO-QUBIT STATE 

In this section, we discuss the optimal time required to reach a two-qubit maximally 
entangled state starting from a pure product state and some partial entangled state. Further, 
we will illustrate how the presence of initial entanglement helps us to minimize the waiting 
time. 

Let the initial state |^E'(0)) of two qubits be 

l^i?) = v^|0)|l) + v/r^|l)|0), (8) 
where < p < 1/2. A general two qubit Hamiltonian can be written as 



3 3 3 

-ff = ^ aiffi ® / + ^ Pjl ® cTj + ^ lijO-i ® cTj, (9) 

i=l j=l i,j=l 

where ai, f3i are real numbers and 7 is a real matrix and cr's are the usual Pauli matri- 
ces. Supplementing the evolution operator with local unitary operators, we can rewrite the 
Hamiltonian as H = Y'jk=i f^kC^k ® cTfc- where yUi > /i2 > /^3 > are the sorted singular 
values of the matrix 7 [5]. Let |^'(0)) evolves under a general two qubit Hamiltonian. So 
the quantum evolution can be written as 



I^b) I^t) = e-'^^/^l'^E) = e-'^^'=i'^^"'=®"'=|^s) = |^+), (10) 

where l'^'^) = :^(|0)|1) + |1)|0)) is a maximally entangled state. The speed with which two 
qubit state evolves in time is given by 

V = = ifi, + (11) 

When the target is a maximally entangled state, then the bound on time is given by 

T > (12) 

2(/ii + /i2)Vl-4Kl-p) 

where Sq = 2cos"^{y/p+ v^l — p)/\^. For the choice of best initial entangled state, we 
have p = po ~ 0.0832 ^]. This state has entanglement E{'^e) ~ 0.413 ebit. To increase 
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the entanglement from 0.413 ebit to one ebit, we must wait for a minimum time given by 

rp _ 0.5911 

~ (m+M2)- 

Next, we ask: does the presence of initial entanglement help to reduce the waiting time? 
The answer is indeed yes. Instead of an initial entangled state, if we start from a product 
state I^E'(O)) = |0)|1) and apply the unitary operator U{T) to reach a maximally entangled 
state (which is our target state), then the bound is given by T > (^^f^^^^ , where 5*0 = 
2 cos~^(| (01|\E'^) I) The minimum time is given by Tmm = • This shows that if we start 

from a product state and want to reach a maximally entangled state, then we will have to 
wait longer. This also illustrates the role of initial condition (or initial entanglement) in the 
brachistocrone problem. 



IV. BOUND ON TIME AVERAGE OF ENTANGLEMENT RATE 

The study of dynamics of entanglement in quantum systems is an ongoing area of research. 
Given a nonlocal Hamiltonian, it is important to know how well one can create entangled 
states. In Ref. [5|, it was found that it is better to start with initial entangled state and best 
initial entanglement is independent of the physical process. In this section, we will show 
that, in a bilocal rotating frame, the entangling rate is in fact a geometric quantity. 

To quantify the entanglement production, one can define the entanglement rate r(t) = 
where E{'^) is some entanglement measure for the state \'^{t)). If we consider von 
Neumann entropy E = — tr(p(t) logp(t)) as our entanglement measure, then we have 

r(t) = -tr(p(t)logp(t)) = -^^^logA„(t) 

n 

N 

= -2 V>'n{t)XUt) logXn{t)KUt), (13) 

n,m=l 

where hnmif) = Im[(a„(t)|(6„(t)|if|am(t))|6m(t))] is called entangling capability of the non- 
local Hamiltonian. In recent times, there have been various attempts to give bound on 
the entanglement rate. For example, upper bound on F in the presence of local ancillas 
was proposed by Childs a/ Wang and Sanders have shown that T < j3 1.9123 for 



any product Hamiltonian H = X ® X ^ where X is a self-inverse Hamiltonian [10|. Ben- 
nett et al have shown that T{H) < cd'^\\H\\, where d is number of Schmidt coefficients, and 
c = OiX) [lH- Childs et a/ have shown that any two product Hamiltonians can simulate each 



other [12j. Bandyopadhyay and Lidar have investigated the effect of noise on the entangling 
capacities of two-qubit Hamiltonians [13;]. Very recently, Lari et al have investigated the 
entanglement rate for two qubits, two qutrits and three qubits using a geometric measure 
of entanglement which holds for multiparticle systems [13] . However, in the literature there 
is no universal bound for entanglement rate. What we plan in this paper is to provide a 
universal bound on the time average of entanglement rate. 

Before, giving a bound on the time average of entanglement rate, we will show that in some 
special cases, the entangling capability is directly related to the speed of the transportation 
of the entangled quantum system. Once we apply the non-local Hamiltonian, the state at a 
later time is given by Eq(l). Suppose, it is possible to perform time-dependent local unitary 
transformations and make the Schmidt basis {|an(^))} and {\bnit))} time independent. This 
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can be imagined as the state of the entangled system in a bilocal rotating frame. Since local 
unitaries cannot change the entanglement, we can define a state 



N 



\^R{t)) = U\t) ® v\t)\^>{t)) = VK{t)K)\hn). (14) 

n=l 

where |a„(t)) = t/(t)|a„) and \hn{i)) = V{i)\bn)- Thus, at any given instant of time the 
entanglement content of \^ nit)) and \^{t)) are the same, i.e., E{^R{t)) = E{^{t)) 

Now, let us look at the (squared) speed at which the state |^ij(i)) evolves in time. It is 
given by 

n ^ ' 

For two-qubits with the general non-local Hamiltonian one has 

^ = 2V^(r^Im[(ao|(6o|/^|ai)|&i)], (16) 
where Ai = p and A2 = (1 — p). Therefore, the speed of the system can be expressed as 

^ pm-p{t)) ^ 4Im[(ao|(6o|i^|ai)|&i)]^. (17) 

Note that h = lm[(ao|(6o|-^|'^i)|&i) and /imax is entangling capability of the non-local Hamil- 
tonian. Therefore, from (17) we have V — 2h. Thus, Knax, the maximum speed of the 
entangled system during quantum evolution is directly related to the entangling capabil- 
ity of two-qubit Hamiltonian. This shows that in a bilocal rotating frame the entangling 
capability is a geometric quantity. 

Next, we will show that one can give a universal bound on the time average of entangle- 
ment for any initial state and general nonlocal Hamiltonians. Consider the evolution of a 
bipartite system |\I/(0)) \^{T)) during an interval [0,T]. We define the time average of 
entanglement f as 

f = ^j\{t)dt = ~Y^j\ogK{t) dK{t). (18) 

Note that during evolution of a composite system the change in the entanglement content 
can never exceed logA^, i.e., 5E = E{T) — E{0) < \ogN, where N is number of Schmidt 
coefficients in the entangled state. This implies that T < ^^^j^. Now, using the geometric 
uncertainty relation for time and energy fluctuation, we flnd 

AO- 

f<21og7V--. (19) 

Thus, the maximal value of time average of entanglement rate depends on the maximum 
Schmidt rank, on the quantum fluctuation in the nonlocal Hamiltonian and inversely on 
the shortest geodesic path. For a given system, the time average of entanglement does 
not depend on the particular Hamiltonian, rather on the fluctuation in the Hamiltonian. 
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There can be infinite number of nonlocal Hamiltonians, yet they all may give same average 
entanglement rate. For two qubits, starting with the initial state as given in (8), the maximal 
time average of entanglement is given by 

r_ = ^^^^v/i-4p(i-p). (20) 

•JO 

It may be mentioned that it is not yet known whether the maximal entanglement and 
disentanglement rates of a nonlocal Hamiltonian are the same, i.e. if T{H) = r{—H)7 We 
would like to point out that even though we do not know complete answer to the above 
question, our result suggests that the maximum attainable value for the time averaged en- 
tanglement rate and time averaged disentanglement rates are identical. The maximum of 
time average entanglement rate is given by Tmax{H) = 2\ogN^^. Since it depends on the 
quantum mechanical uncertainty, we have AH = A{—H) and because of the law of reci- 
procity of quantum mechanical probability (5*0 is same during entangling and disentangling 
evolutions) we have Tmax{H) = Tmax{—H). This is another non-trivial consequence of our 
geometric approach to the problem. 



V. MINIMUM TIME FOR COMPOSITE HAMILTONIAN 

Now, we ask what happens to the minimal time if we consider a time evolution under one 
Hamiltonian instead of another. In particular, if the later is build up from a composition of 
Hamiltonians then how does the minimum time behave. The answer to this question may 
be useful in quantum computing algorithms, where one tries to design sequence of unitary 
operators via the application of suitable Hamiltonians. It is known that if we have the 
ability to perform evolution under a Hamiltonian H (say) and have the ability to apply 
unitaries Uk, it is possible to simulate evolution according to a composite Hamiltonian like 
H' = Y2f^ akUhHUl, where a^s are real numbers 1^. Under such an evolution we will have 



r>^cos-M(^(0)|t/(T)lM/(0))|, (21) 

where U{T) = exp{—iH'T / h) and AH' is the quantum fluctuation in the composite Hamil- 
tonian. Suppose, that the initial and the final states are fixed, i.e., we start with same 
initial state and want to reach a desired target state with the new Hamiltonian H'. Then, 
the shortest path 5*0 is again the same. Now, using the convexity property of quantum 



mechanical uncertainty [15|, we have 



AH' = A{Y,o^kUkHUl) < J2(^kA{U,HUl). (22) 

k k 

This tells us that mixing different Hamiltonians decreases the quantum mechanical uncer- 
tainty. If we define Hk = UkHUl, then the minimum time for the new Hamiltonian satisfies 

- 2AH' - 2j:,akAHk ^ ^ 

This also gives a composition law for the minimum time for entangled brachistocrone 
problem under composite Hamiltonian. To illustrate this, consider a composite Hamiltonian 
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H' that consists of two Hamiltonians Hi and H2. If wc evolve the system only under Hi 
then the minimum time to reach the target entangled state will be Ti^^^ = ^^^^ . Similarly, 
if we evolve the system only under H2 then the minimum time to reach the target entangled 
state will be T2^^^ = j^^- Now, if we evolve the system under H' — aiHi + 0:2 -ff2 then the 
minimum time to reach the target entangled state will be 

T T 

rp -'- mini -'- min2 {0A\ 

-'-min — rp I ^ rp ■ V-^^J 

Otl-Lmin2 + Ot2-Lmin\ 

To illustrate this composition law consider the initial state as given in (8). Now suppose, 
we apply the Hamiltonian Hi = /xicri (Ji and evolve the system for sometime to reach 
the target entangled state. In this case, the minimum time we need to wait is given by 
Ti . — If we apply instead H2 — 112(^2 ® ^2 and evolve the system to reach 

the target entangled state, then the minimum time we need to wait is given by T2.^^^.^^^ = 
Now, we apply H = aiHi + Q;2-f^2 and evolve the system for sometime to 

2/t2Vl-4p(l-p) 

reach the same target entangled state. In this case, the minimum time we need to wait is 
given by T^jn = , ■ This later quantity indeed satisfies the composition 

law as given in (24). In fact, with a suitable controlled experiment, one can even test the 
composition law for the minimum time. 



VI. CONCLUSION 

We have investigated the brachistochrone problem for entangled states. This answers 
the question that we have raised in the beginning, namely, given an arbitrary product 
state and general non-local Hamiltonian, what is the minimum time for creating a desired 
entangled state? We have shown that the minimum time depends inversely on the speed of 
the evolution of the quantum system. This is independent of the particular entanglement 
measure and holds even for multiparticle entangled systems. We have shown that for two 
qubits the presence of initial entanglement does help in minimising the waiting time. In 
particular, we find that if we are given an option of starting from a product state and an 
entangled state, and wants to reach the maximally entangled state, then it is better to start 
from some already initial entangled state. Because, in the later case the time required is 
less. We have shown that the time average of entanglement rate depends on the Schmidt 
rank and on the fluctuation in the non-local Hamiltonian. We have also shown that in a 
bilocal rotating frame the entangling capability is also directly related to the speed of the 
evolution of the entangled states. This shows that the entangling capabihty of a non-local 
Hamiltonian is actually a geometric quantity. Furthermore, we have proved a composition 
law for the minimum time when the system evolves under composite Hamiltonians. We hope 
that these findings will be of interest in understanding the entanglement rate and entangling 
capabilities of non-local Hamiltonians and in real experimental situations. In future, it will 
be interesting to see if all these findings also hold for qudits and multiparticle systems. 
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